Quantum simulator for the Hubbard model with long-range Coulomb interactions 

using surface acoustic waves 
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A practical experimental scheme for a quantum simulator of strongly correlated electrons is pro- 
posed. Our scheme employs electrons confined in a two dimensional electron gas in a GaAs/AlGaAs 
heterojunction. Two surface acoustic waves are then induced in the GaAs substrate, which create 
a two dimensional "egg-carton" potential. The dynamics of the electrons in this potential is de- 
scribed by a Hubbard model with long-range Coulomb interactions. The state of the electrons in 
this system can be probed via its conductance and noise properties. This allows the identification 
of a metallic or insulating state. Numerical estimates for the parameters appearing in the effective 
Hubbard model are calculated using the proposed experimental system. These calculations suggest 
that observations of quantum phase transition phenomena of the electrons in the potential array 
are within experimental reach. 

PACS numbers: 03.67.Lx, 71.10.Fd, 74.25.-q 



Since the pioneering experiments of Greiner et ai, 
where bosonic atoms were trapped in optical lattices to 
create a Bose-Hubbard model (jy], quantum simulators 
have attracted a large amount of attention. The interest 
is twofold. The first is that such devices are hoped to offer 
an alternative method for studying quantum many-body 
systems, which appear in many branches of physics, yet 
in many cases no efficient and reliable means are available 
to obtain quantitative information about them. The sec- 
ond is from the perspective of quantum computing, where 
the technology developed for such devices are hoped to be 
useful for realizing a scalable quantum computer. Indeed, 
many proposals for using optical lattices in quantum in- 
formation processing devices exist 0. 

In this paper, we propose another type of quantum 
simulator that may be used for studying fermionic par- 
ticles interacting via a Hubbard-type interaction. The 
Fermi-Hubbard model, as opposed to the Bose-Hubbard 
model, is particularly interesting due to its significance in 
condensed matter theory, ranging from its original incep- 
tion in understanding metal-insulator transitions, to its 
recent reincarnation in high-tcmperature superconduc- 
tivity. In particular, it is still controversial in whether 
the Hubbard model supports a c?-wave superconducting 
phase ■ It is also the basis for gate operations in elec- 
tron spin-based quantum computing [J. Recently there 
have been many advances in the trapping of fermionic 
atoms in optical lattices |5|| . One of the features of the op- 
tical lattice approach is that the effective form of the in- 
teractions is a contact interaction that is extremely short 
ranged |^. Since in a real system displaying Hubbard 
dynamics the interactions are of a long range form origi- 
nating from the Coulomb interaction, such effects may be 
difficult to be taken into account in an optical lattice ap- 
proach. One may speculate that to simulate a Coulomb 
interaction, a real Coulomb interaction is necessary in 




FIG. 1: (Color online) A sketch of the proposed experimental 
structure. Two IDTs launch SAWs (green) at right angles to 
each other forming a moving two-dimensional "egg-carton" 
potential. Light blue regions indicate the regions where the 
2DEG is present. Part of the 2DEG channel is shallow etched 
(yellow), to partially deplete the electron density. Source (S) 
and drain (D) ohmic contacts (dark blue) are formed at the 
ends of the 2DEG. 



the system carrying out the simulation. 

Our proposed experimental setup is shown in Fig. 
n We start with a standard modulation-doped 
GaAs/AlGaAs heterojunction, which forms a two dimen- 
sional electron gas (2DEG). The 2DEG is formed into a 
channel as shown in the figure, with source and drain 
ohmic contacts placed at the two ends. Between the 
source and drain contacts, a region of low electron density 
is formed via locally raising the conduction band, which 
may be done by a shallow etching procedure. The ohmic 
contacts thus probe the conductance properties of this 
region. Outside of the 2DEG mesa, interdigitated trans- 
ducers (IDTs) are placed at right angles to each other. 
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Surface acoustic waves (SAWs) are launched by apply- 
ing a high frequency AC voltage to the IDTs, forming an 
interference region at the center of the device. Due to 
the piezoelectric property of GaAs, an electric potential 
is induced following the shape of the SAW modulation. 
The net potential due to the SAWs in the central region 
is then equal to the sum of the individual SAW poten- 
tials. Assuming both IDTs produce identical SAWs, the 
total potential at a particular time is then 



V{x, y) = Vb (cos(fca:;) + cos(fcy)) , 



(1) 



where Vb is the amplitude generated by the SAWs, k = 
AW is the wavenumber of the SAWs, with ^saw 
being the SAW velocity (usaw ~ 2860 m/s in GaAs), and 
/ the SAW frequency. This is an "egg-carton" shaped 
potential that moves in the diagonal direction towards 
the drain electrode. We assume that the approximate 
size of the interfering region that can be produced is 50 
/imx50 /xm. The depth of the potential can be controlled 
by varying the amplitude Vb of the SAWs. 

Such a moving potential may be used to transport elec- 
trons from the source to drain electrodes in a similar way 
to that observed in Ref. 0- As the SAW moves, it car- 
ries the trapped electrons with it in the local minima of 
the SAW (the acoustoelectrical current). If we assume 
that n electrons are trapped in each minima, the acous- 
toelectrical current is then given by / — nef , where e 
is the electron charge. A crucial difference between the 
current setup and the experiment of Ref. 7] is that the 
length of the region between the source and drain elec- 
trodes was only of the order of a single SAW wavelength 
in Ref. 0- Therefore, only one or two potential traps 
are realized in the central region in their case. Assum- 
ing a SAW frequency of / « 20 GHz in both directions 
(corresponding to a wavelength of A = 0.14 /im), we have 
^ 10^ potential minima in the central region. The effec- 
tive lattice size will be governed by the coherence length 
of the sample. Long coherence lengths are now achiev- 
able due to advances in modern fabrication techniques 
such that large sections of the lattice will interact coher- 
ently. For example, mobilities ^ corresponding to mean 
free paths of Z„,ip — vp^im* je « 11 /im at T=0.1 K and 
electron densities of rid = 10^*^ cm~^ were reported in 
Ref. 0, where ?;f is the Fermi velocity (wi? w 2 x 10^ 
cm/s in GaAs), and to* is the electron effective mass 
(to* — 0.067TOe in GaAs). The properties of the current 
flowing from the source to the drain are therefore de- 
pendent on the collective effect of the lattice of electrons 
formed by the SAWs. 

To model the electrons in this interference region, let 
us start with the general Hamiltonian 
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where V (z) is a triangular potential well associated with 
the heterojunction, and ipiri^) ^re the fermion field op- 
erators for electrons of spin a. The Coulomb potential 
Uc{x, x') experienced by the electrons in the 2DEG will 
be strongly screened due to the high density of electrons 
in the (5-doping layer located near the 2DEG. These elec- 
trons are separated from the heterojunction by a distance 
d, which may be controlled by the thickness of the spacer 
layer. We therefore take the form of the Coulomb poten- 
tial to be 



+ ^E/ d^xSx'i^l,{x')^l{x)Uc{x,x')^,{x)i^,,{x'), 



Uc{x,x') = 



1 



Ane \x — x'\ 



(2) 
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\x — X 



^J{x - x'Y + iy- y')'^ + [z + z' + 2dy 



where e is the permittivity (e « 13eo in GaAs). Equation 
was derived assuming a metal plate in the x-y plane 
a distance d from the electrons. 

The single particle Bloch wavefunction solutions of TL 
are simply Mathieu functions, from which we may con- 
struct a Wannier basis ^ 0| . There is a localized Wan- 
nier function w{x — Xj) for every local minima in the x-y 
plane, labeled by the index j. Making the transformation 
Tpaix) = Cja-w{x ~ Xj), whcrc CjCT is a fermion annihi- 
lation operator associated with the site j, we obtain the 
Hamiltonian 



CT ij 



EE^ 



CkaCla' 



(3) 



ctct' ijkl 



where n^o- = cl^Cia and 



tij = I d^xw*{x — Xi) 



2m 



- + F(x, y) w(x — Xj) 
(4) 



Uijki — j d^xd^x'w*{x' — Xi)w*{x — Xj)Uc{x,x') 

xw{x — Xk)w{x' — xi). (5) 

If one chooses to discard all but the largest matrix el- 
ements in Eqs. Q and one obtains the standard 
Hubbard model (in standard notation, t = —tij, for i,j 
nearest neighbors and U = Uuu ). It is an important fea- 
ture here that the full effect of the Coulomb interaction 
is present in the Hamiltonian including long-range 
interactions and correlated hopping terms. These are 
difhcult to take into account from a calculational point 
of view, but are naturally included in the current setting. 

Let us evaluate the integrals Eqs. and (0) for 
the proposed experimental parameters. The results are 
shown in Figs. [21 and |3| First, we find that the only 
hopping directions that are allowed are those along the 



3 





1 


1 


1 














toy" 








"8. , 
















- 
















" " " ' ■ - - . . . _ 

1 


1 






i ------- 



0.3 

V„(racV) 



FIG. 2: Various Hubbard parameters as a function of the 
SAW potential height, t is the nearest neighbor hopping, t' is 
the next-nearest neighbor hopping, V — Uijji is the nearest 
neighbor Coulomb integral, C2 = Unjj is a pair correlated 
hopping term, with i,j nearest neighbor sites (see inset of 
Fig. 1^. The lattice sites involved in the various terms are 
shown in the inset of Fig. El The chosen parameters are a 
SAW frequency of / = 20 GHz and a screening distance of 
d = 10 nm. Effective mass and permittivity values were taken 
to be those in GaAs. Error bars are due to the Monte Carlo 
integration procedure used. Lines are to guide the eye. 



lattice axis directions, i.e. all diagonal hoppings are zero. 
Therefore, the next-nearest neighbor hopping term t' is 
a two lattice site hopping term as shown in the inset 
of Fig. 13 This is simply a result of the Wannier ba- 
sis that was chosen to write the Hamiltonian Eq. 
and is due to the symmetry of the square lattice that 
is currently being used. All hopping terms decay with 
a roughly exponential dependence on the SAW potential 
height Vo . The on-site Coulomb term U increases with Vq 
as expected due to the greater localization of the Wan- 
nier functions for larger potentials. The nearest-neighbor 
Coulomb term V (i.e. the coefficient of the term rii^nj^i 



with i,j nearest neighbors) 
value of V{Vo 00) 



approaches the asymptotic 
1 

' We 



. i/v/TTIrf^A 

8 /ieV, consistent with the numerical calculation, 
also calculate one of the correlated hopping terms de- 
noted as C2 in Fig. 2 which is the coefficient of the term 
c'j^c'jf^Ci-^Cii where i,j are nearest neighbor sites. Due 
to the small amplitudes many of the correlated hopping 
terms are numerically difficult to calculate. This corre- 
lated hopping term appears to have the largest amplitude 
of these. 

One expects that due to the Hubbard nature of the 
electron interactions in Eq. Q that as the SAW poten- 
tial is raised, at some point a metal-Mott insulator tran- 
sition will be reached [l^. Figure^ shows the trajectory 
of the Hubbard parameters in the space (U/t,V/t,t' /t) 
as the SAW potential is increased. On this trajectory we 
have superimposed a magnetic phase transition bound- 
ary in this parameter space as obtained from mean field 
theory, following Refs. [ill fl^ . We use the Hamilto- 
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FIG. 3: The on-site Coulomb integral as a function of the 
SAW potential height. The inset shows the connecting points 
for the various Hubbard parameters calculated in Fig.|21 Lines 
are to guide the eye. 




FIG. 4: The trajectory in the space {U/t, V/t, t' /t) as the 
SAW potential is increased. The phase boundary between 
the paramagnetic metal (PM) and antiferromagnetic (AF) re- 
gions, as well as the boundary between the antiferromagnetic 
and ferromagnetic (F) regions (both insulating) as obtained 
from mean field theory at filling n = 1.3 is also shown. The 
stars indicate the points where the parameters pass through 
the boundaries. Dashed lines show the projection of the tra- 
jectory onto the planes spanned by the {U /t,V /t,t' /t) axes. 



nian Eq. ||2Jl, with all parameters other than U,V,t' and 
t set to zero. Self-consistent solutions are found for para- 
magnetic, antiferromagnetic, ferromagnetic, and charge 
density wave states, and the phase diagram is mapped 
by taking the lowest energy solution. We see that the 
system crosses over from a paramagnetic metal phase 
to an antiferromagnetic phase, then to a ferromagnetic 
phase as the SAW potential is increased. We find that 
the effect of the two-site hopping t' is qualitatively simi- 
lar to the more conventional diagonal hopping term, with 
good agreement between our results and those obtained 
by a diagonal coupling given in Ref . • Although it is 
unclear whether the paramagnetic to antiferromagnetic 
transition as calculated in mean field theory corresponds 
precisely to the Mott transition point, our results sug- 
gest that our parameter values are in the correct range 
to observe quantum phase transition phenomena. 
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This transition will be reflected in a modified band 
structure in the central region of the setup shown in 
Fig. nj In Fig. we show the band structures due 
to two SAW intensities corresponding to the metallic 
and Mott insulating phases. First consider the weak 
SAW modulation case (Fig. The periodic poten- 

tial Eq. 1^ produces a conduction band as shown in the 
figure. Due to the weak modulation, the band struc- 
ture corresponds to a metallic state, i.e. there is no gap 
between the ground states of successive electron num- 
bers. The band structure due to the strong SAW mod- 
ulation is shown in Fig. [SJs. Ground states with succes- 
sive electron numbers N = separated from 
each other by a charge gap A. Now consider applying a 
small source-drain voltage Vds- In general this will add 
a dc-current /dc component to the acoustoelectrical cur- 
rent. This dc-current probes the excitations or the phase 
of the Hubbard model. Considering the Mott insulator 
regime (Fig.Eb), when a small bias voltage Vds such that 
eVds, ksT < A is applied, the dc-current is blocked as a 
consequence of strong correlations. On the other hand, 
in the metallic phase no gap opens and Ohm's law should 
hold /dc oc Vds (Fig-Et)- The phase transition is reached 
by tuning the 2DEG electron density n = N/A where A 
is the area of the SAW lattice. The density can be con- 
trolled by a backgate voltage and/or the SAW potential 
strength Vq. 

The absence or existence of fluctuations in the oc- 
cupation numbers fii can be probed by low-frequency 
shot noise. The shot noise due to /dc is supressed in 
a macroscopic conductor However, in an insulat- 

ing phase the fluctuations in the SAW driven acous- 
toelectrical current should still be observed due to the 
quantized bunches of electrons per SAW-potential min- 
ima. Assuming zero, one, or two electrons per site, the 
shot noise is proportional to f{{5nif') with {{SfiiY) = 
{ni){l - {fii)) + 2{ni^nii) where fii = fii-^ + fni 
The complete absence of doubly occupied and empty sites 
is only predicted at half-filling in the limit of infinite U 
[lij . Therefore, at half-filling and deep into the Mott 
insulator regime the insulating phase is characterized by 
the disappearance of shot noise. Off half-filling, the elec- 
tron number on a given site is free to fluctuate and we 
expect a finite shot noise. Noise measurements in combi- 
nation with conduction measurements therefore can clar- 
ify the role of fluctuations in the occupation numbers in 
the insulating regime. 

In order to observe a quantum phase transition, we 
require fc^T to be smaller than the typical energy scale 
of the Hamiltonian, which is set by the hopping inte- 
gral t 19] . We expect that a phase transition will occur 
for relatively small values of the SAW potential, where 
t K, 0.05 meV. This corresponds to temperatures less 
than T « 0.6 K, which is an experimentally accessible 
temperature. Another consideration is due to the fact 
that we are using a moving potential which remains in 
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FIG. 5: (Color online) Band structure (cut along the diagonal 
in Fig.0of the Hubbard simulator coupled to source (S) and 
drain (D) reservoirs. The metallic phase is depicted in a). 
The energy to add or remove electrons via the source and 
drain electrodes is gapless which means that the Fermi level 
/i (dashed line) of the reservoirs lies within the conduction 
band (blue). The Mott-insulator phase is shown in b) for 
filling n — 1. Changing the number of electrons in the lattice 
region is associated with a charge gap A which suppresses 
dc-transport if fcsT, eVds < A. This is indicated by a lifted 
conduction band (red) in the bulk of the lattice. 



the central region of Fig. ^ for a time of approximately 
L/vsA-w, where L is the length of the central region. This 
time corresponds to approximately 15 ns for L 50 /im, 
which must be larger than the tunneling time for the elec- 
trons in the lattice. As the tunneling time is of the order 
h/t 10 ps, there is time for many tunneling events and 
therefore equilibration should not be problematic in this 
case. 

A simple variation on the current proposal is to inde- 
pendently vary the SAW potential heights in the x and 
y directions, which changes the hopping amplitudes 
in these directions. By making the SAW potential much 
larger in one direction, we may smoothly evolve the sys- 
tem from a 2D lattice to a ID chain. Another possible 
variation is to add a random potential R{x, y, t) to the 
potential Eq. This is easily implemented by apply- 
ing a random amplitude modulation voltage to the IDTs. 
The effect of such a term is to randomly distribute both 
the chemical potential /i and hoppings tij from site to site 
in the Hamiltonian Eq. ||3Jl. This is a Hubbard- Anderson 
model and has been considered in studies of interacting 
disordered systems |20j . The random potential can be 
characterized by two parameters, the standard deviation 
of the amplitude which controls the degree of disorder, 
and the correlation J R{x -I- Ax, y -I- Ay, t)R{x, y, t)dxdy 
which controls the density of the "impurities" . By si- 
multaneously exciting higher harmonics of the IDT, and 
length modulating the amplitudes of these harmonics 
randomly, the correlations of R{x,y,t) may be varied. 
In this way we may study metal-insulator transitions 
brought on by the effects of disorder, as well as corre- 
lations. 
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